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We consider the entropy associated with the large-scale structure of the Universe in the linear 
regime, where the Universe can be described by a perturbed Friedmann-Lemaitre spacetime. In 
particular, we compare two different definitions proposed in the literature for the entropy using a 
spatial averaging prescription. For one definition, the entropy of the large-scale structure for a given 
comoving volume always grows with time, both for a CDM and a ACDM model. In particular, 
while it diverges for a CDM model, it saturates to a constant value in the presence of a cosmological 
constant. The use of a light-cone averaging prescription in the context of the evaluation of the 
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I. INTRODUCTION 

In the standard cosmological approach, the Universe is 
described by a homegeneous and isotropic solution of Ein¬ 
stein field equations, known as Friedmann-Lemaitre (FL) 
spacetimes [T]. These solutions are expected to describe 
the Universe smoothed on cosmological scales. While this 
spacetime is easily identified in the early universe, matter 
clusters and structures grow under the effect of gravity so 
that the distribution of matter in our late time universe 
exhibit large inhomogeneities. 

Over the past years significant activity has been de¬ 
voted to the definition of averaging procedures [MS] in 
order to construct a notion of a coarse-grained spacetime. 
Irrespective of the question of whether this procedure 
could explain the recent acceleration of the Universe, 
averaging methods are interesting since let us compare 
the evolution of a system described at different scales. 
In such a coarse-graining, information about the micro¬ 
scopic behaviour of the system is lost, which is at the 
origin of the notion of entropy. In essence, entropy esti¬ 
mates the number of micro-states that correspond, after 
averaging, to a given macro-state. 
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The definition of gravitational entropy is still an open 
debate. While a suitable definition has been given in 
the context of the thermodynamics of stationary black 
holes [6] , a well motivated and universally accepted ana¬ 
logue has yet to be found in the cosmological context. 
With the evolution of the Universe, structure grows and 
the Universe becomes more and more inhomogeneous. In 
order for the second law of thermodynamics to hold, the 
gravitational field itself shall carry entropy. It was ar¬ 
gued in I3i that it has to be defined from the free grav¬ 
itational field and thus be related to the Weyl tensor; 
more recently it was extended to a definition [S] based 
on the Bel-Robinson tensor. This latter proposal reduces 
to the Bekenstein-Hawking entropy when integrated over 
the interior of a Schwarzschild black hole and increases 
as inhomogeneities grow. 

While the concept of entropy arose from equilibrium 
thermodynamics, it has been also thought of as a mea¬ 
sure of information. In terms of information theory the 
Kullback-Leibler divergence m, for two probability dis¬ 
tribution functions p and q, is defined by 

DKL(pk) = ^In^^ = J p{x)\n^^dx , 

and quantifies the amount of information lost when the 
data (p) is represented by the model (q). In cosmology, 
it was used in order to decide whether two cosmological 
models can be distinguished given a set of observational 
data [TT] . 
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Since the cosmological model of structure formation 
predicts the distribution of the density field (as a random 
variable), it has been proposed [12], in the context of 
averaging, to adapt the Kullback-Leibler divergence to a 
definition of relative information entropy that quantifies 
how the actual density field p is different from its spatial 
average (p)® on a spatial domain V of proper volume 1^, 
namely ^ 


^Rl; 


V 


1 


V-D Mpi 


pin 


(P) 


V 


( 1 ) 


V 


with Mpf = SttG. It was conjectured that this function 
is increasing with cosmic time, which was checked for 
linear perturbations of a spatially Euclidean FL space- 
time for a CDM model m and for the comparison of 
a Lemaitre-Tolman-Bondi (LTB) spacetime to its aver¬ 
age [T^17j . This is indeed a key property to consider 
Eq. Q as a valid definition of entropy. 

This definition lies solely on the density field. In princi¬ 
ple, we want to compare a spacetime (Ai , g) to its average 
(M, g), and evaluate the quantity of information that has 
been lost. The previous definition assumes that Ai can 
be foliated by a family of spatial (or null) hypersurfaces, 
which is related to the choice of the averaging procedure. 
It means that (i) the average depends on the choice of 
the slicing. The expression in Eq. Q involves an integral 
over quantities in the two spacetimes. It means that it 
implicitly involves a mapping between A4 and AA so that 
(ii) it may have a gauge dependence, (hi) In principle, 
it also depends on the averaging procedure, whether it is 
spacelike [2] or along the light-cone |5]. The expression 
0 also depends solely on the matter distribution. After 
averaging, the matter distribution is indeed homogeneous 
but it may not be isotropic, as e.g. shown in Ref. m, 
so that (iv) it may not capture the fact that the back¬ 
ground spacetime may not be FL. Indeed, the geometry 
of the spacetime Ai is characterized by its Riemann ten¬ 
sor R^vpcr, that can be split as a Ricci contribution, 
and a Weyl part (that can further be decomposed 

as an electric and magnetic parts, and B^^). The 
Ricci part is constrained by the matter distribution, via 
the Einstein field equations, but the knowledge of the 
density field alone does not allow one to reconstruct the 
Weyl part. As an example, consider AI as a perturbed 
FL universe; it has non-vanishing Rp^, Ep^ and Bp^. Af¬ 
ter averaging (p)x) is homogeneous. If TM is a FL universe 
then only Rp,^ is non-zero while if AI is a Bianchi I uni¬ 
verse both Rpy and Epi, are non-zero, while they have the 
same (p)!?. It means that one part of the difference be¬ 
tween the spacetimes is not included in the definition Q . 
Part of this information is contained in the two scalars 

for definition) 


lie 


Cp,paC>^''p'^ and C*p,p,CP''P^ (see Sec.^ 
constructed from the Weyl tensor, which is indeed at the 


heart of the proposals [^. While not transparent in 
the definition Q, it can actually be shown (see Ref. [T3|) 
that, for perturbations around an FL spacetime, this for¬ 
mula reduces in part to some combination of the Weyl 
scalars. 

The paper is organized as follows. In Sec.jnjwe define 
the foliation of our spacetime, the two averaging proce¬ 
dure we shall consider and introduce linear perturbation 
theory. In Sec. |III| we present the different definitions of 
the gravitational entropy used, while in Sec. |IV| we inves¬ 
tigate their time evolution. Finally, the results are dis¬ 
cussed and compared in Sec.jVj Appendix [Ajpresents the 
dynamics of the background spacetime, and Appendix |B| 
summarizes the definition of gauge invariant degrees of 
freedom. 


II. AVERAGING PROGEDURES 


Averaging procedures rely on a choice of ob¬ 
servers/foliation of spacetime, described in Sec. IIA In 


the following, we consider two procedures base d either 
on spatial sections or on null sections (Sec. IIB). 


A. Spacetime foliation 


To define our formalism let us introduce a 1-1-3 split¬ 
ting of the Universe [TS] associated with a general ref¬ 
erence timelike congruence that defines a class of 
observers and the relative foliation of spacetime. The 
3-dimensional spacelike hypersurfaces normal to can 
then be defined by the equation S{x,t) — Sq = 0, with 
S{x,t) a scalar field and Sq a constant. Then 


nu. = - 


dpS 


{-dpSd.s gp-y/^^ 


( 2 ) 


is normalized as Upn^ = —1. This allows us to define 
hp^, the projector on these hypersurfaces, as 

hpv = g^v + npU^ (3) 

which satisfies by construction hpphP = hpi, and hp^nP = 
0. Furthermore, one can define the expansion 0, shear 
and vorticity ujp^ of the flow as 


0 ,,. = hT,h?y, 
.0 

They are explicitly given by 


- 


•' fJ.l' 


+ W 




0 = ^puP, 


(4) 

(5) 

( 6 ) 






(7) 


^[iiy — [ot'^0] ■ 


^ We add a factor IjMpi with respect to the definition proposed 
in m to obtain a dimensionless relative information entropy. 


( 8 ) 
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Indeed, the assumption of Eq. ([^ implies that the vor- 
ticity strictly vanishes, = 0. 

In practice, perturbations grow significantly only dur¬ 
ing the matter-dominated era, so that one can restrict 
the analysis to dust-filled universes, eventually with a 
cosmological constant. In such a situation one can pick 
up a foliation dehned by a congruence corresponding 
to the four-velocity of a geodesic observer, which acci¬ 
dentally coincides with the four-velocity of comoving 
observers, i.e. 

The shear tensor can be then expressed as 


where M is the 4-dimensional spacetime. Vs is the vol¬ 
ume of the 2-sphere embedded in the light-cone, defined 
by the requirement that (l)vb,So = 1, and g is the deter¬ 
minant of the four dimensional metric 


C. Linear peturbation theory 

The standard scalar-vector-tensor decomposition [T] of 
a perturbed FL spacetime have metric components 


- ^ UL' r, : 


and it follows that the scalar shear takes the form 

= l( 0 (; 0 " - ^02 

r\ y jjL V (J, q 


(9) 

o o 

(10) 



= —2a^a 


= a 




,(13) 


B. Spatial and light-cone averaging 

On one hand, we consider a spatial averaging proce¬ 
dure [2] , entirely based on a slicing of spacetime by spatial 
hypersurfaces. The spatial average of any scalar quantity 
A on a domain T) is then defined as 

{A{r],x))T:, = ^ ( 11 ) 

where Vx> is the volume of the domain, defined by the 
requirement that (l)x> = 1, and h is the determinant of 
the induced metric on the averaging hypersurface. 
Such a spatial average is associated with the general ref¬ 
erence timelike congruence of Eq. if and only if 
the average is performed in the gauge where S{ri,x) is 
homogeneous (see Refs. [2f)U22] l. 

On the other hand, cosmological observations are usu¬ 
ally restricted on the past light-cone, since most of the 
relevant signals are of electromagnetic origin. Hence, 
when we look to cosmological observables, the averaging 
procedure should be possibly referred to a null hyper¬ 
surface coinciding with our past light-cone or to the null 
surface obtained from the intersection of our past light- 
cone with some fixed-time spacelike hypersurface. Let 
us consider this latter possibility. Following Ref. [5] we 
obtain that the averaging of any scalar A{r], x) over the 
2-sphere embedded in our past light-cone, defined by a 
null scalar V{r],x) (i.e., such that dfjVd^V = 0) equal 
to a constant, and corresponding to its intersection with 
the spacelike hypersurface S{ri,x) = Sq, is given by 

{A{g, x))vo,So = ^ [ V^S{Vo- V)S{S - So) 

Ks Jm 

A{g,x)\df,Vd^S\d'^x, (12) 


with Dij = didj — ^SijA. We then have 4 scalar degrees of 
freedom {a, /3, ip and E), 2 transverse vectors {Bi and Xi 
with d^'Bi — 0, d^Xi = 0) with 4 degrees of freedom, and a 
traceless and transverse tensor (hij with d^hij = 0 = h]) 
with 2 degrees of freedom.^ 

Let us stress that first order perturbation theory is suf¬ 
ficient to obtain the general expression for the shear in 
Eq. ( [To| ) up to second order. Since second order pertur¬ 
bations contribute only to third or fourth order to cr^ (see 
Appendix [b|). 

In the following we shall use the synchronous gauge 
and neglect vector and tensor perturbations. We then 
have 


ds^ = {— d?7^ -I- [(1 — 2ip)Sij + DijE] dxMx-^ } . (14) 


It is clear from Eqs. ( B1||B3 ) that we can then write the 
Bardeen potentials 4/ and 4> as 


1 TV TV W'' 

d! = ip+-AE+—E', $ = -—E'-, (15) 

^6 2 ’ 2 2 ' ^ ’ 

where the prime denotes the derivative with respect to 
conformal time and % — a'/a. 

Let us now introduce the Weyl tensor defined as 

CfiuXp — -j- — {gppRux T gvxBpp g^xBup 

9i^pRpx) d~ g iSpxgi'p 9pp9vx) R 1 (16) 

where Rp^Xp and Rpv are the Riemann and Ricci ten¬ 
sors, while R is the Ricci scalar. We can then define the 
dual of the Weyl tensor as where 

gapr-y = xf-^f-apry IS the four dimensional volume ele¬ 
ment. 


^ In general Up, which defines a general reference flow, and Up, 
which defines the four-velocity of the observers comoving with 
the matter, may be different (see Ref. m for details). 
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The definition of gauge invariant degrees of freedom are summa¬ 
rized in Appendix [B| 
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In terms of the Bardeen potentials, in any gauge 
and for vanishing anisotropic stress, the Weyl scalar 
takes the simple form 


^ A,$ . (17) 

Let us relate the contraction of the Weyl tensor in 
Eq. ( [l7| to the shear of a free-falling observer given in 
the synchronous gauge by 

a" = ^D,,E'D^^E ', (18) 


which shows that cr^ and Cpi^xpC^’'^^ are not indepen¬ 
dent quantities. If we consider a general ACDM model 
we have 


= Q.^2^ V^$(?7,f) -k ^4>(7?„,g), (19) 

E{ri,x) = ~ 3 -^ 2 ^ (20) 


Using Eq. (201 and the background dynamics (see Ap¬ 
pendix A), we obtain the following relation which con¬ 
nects the shear to the gravitational potential 


a 


2 


9 \ao'H'^flmo 


2 

+ Aj$'] 


X [-HA^^-h A^$'] . 


( 21 ) 


The relation between the Weyl tensor in Eq. ( [T7| ) and the 
shear of a free-falling observer is then given by 


Cp,xpC^’'^P 


8 2 / 1 
9 \aoHlilmo 



1 fi 8 

-AAA^$' - —^ . (22) 

a^H ^ ^ ^ ’ 


In a CDM model, the gravitational potential is constant, 
so that the contraction of the Weyl tensor and the shear 
are related simply by a time dependent factor. In a 
ACDM model, one needs to include the terms arising 
from the decay of the gravitational potential. 


III. DEFINITIONS OF THE ENTROPY 
A. Definition from the density field 

As discussed in the introduction, the first idea to define 
a relative entropy between two spacetimes m followed 
the definition of the Kullback-Leibler divergence m in 
information theory. It allows one to quantify whether 
two density fields p and {p)t) can be distinguished and is 
defined^ by Eq. Q. 


In the following we consider the spatial averaging prescription 
< .... >x) in our definition, we will show in Sec. IV that this is 
indeed the right averaging prescription for the evaluation of the 
entropy. 


When working in perturbation the density field can be 
expanded as p = +P^‘^\ so that at leading order 


5ri,-d 

Ud 


1 ((pW)"A-(p«) 


2Mpi 


P 


( 0 ) 


(23) 


As can be seen from this equation, the relative infor¬ 
mation entropy can be obtained at leading order using 
only first order perturbation theory. In particular, this 
is given by the variance of the energy density. 

Eq. (23) is valid independently of the matter content of 


the Universe. Following our introductory considerations, 
the density field that constraints the Ricci part of the 
Riemann tensor is the total energy density. Therefore, to 
define the entropy we shall use the total density energy 
of the Universe. 


Let us now evaluate Eq. (23) in a ACDM universe. 
In the synchronous gauge, the first order perturbation 
of the energy density is given by the Poisson equation 
(this corresponds to the matter perturbation because the 
cosmological constant cannot be perturbed by definition) 


As a consequence Eq. (p3| becomes 


^RR-D _ 2 Mpi 


Vx> 


L 


((V2ci>)%-(V2$)| 


(24) 


(25) 


Let us now use the expression of Eq. HH) for the Weyl 
tensor in a perturbed FL metric to rewrite Eq. (25) in a 


useful form. After some simple algebraic manipulations 
we obtain that 


— - j""' 


18-^2 


{Cp^xpC^'^^np 






As shown in Ref. 


(26) 


the term in the second line of 
Eq. (|26|) is related to the so-called kinematical backreac- 


tion Qx, (see Ref. i), given by 


Qx> = g ((0^)x> - (0)1)) - 2(cr^)-D , (27) 

in the CDM case. Therefore, for a free falling observer 
in the synchronous gauge and considering a CDM model 
we can rewrite Eq. (26) in the following way [13] 


^ = ^-Mpi 

Vv 4 


18-^2 


(C'^.Ap 


(28) 


The entropy (26) is the average of a combination of 


scalar quantities and it is gauge invariant under a gauge 
transformation (see Ref. [55] for the possible gauge de¬ 
pendence coming from the averaging prescription). In¬ 
deed, this scalar combination is zero at zero and first 
orders (see Eq. (23l), and therefore gauge invariant at 


leading order under a gauge transformation. 
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B. Definition from the Beii-Robinson tensor IV. TIME EVOLUTION OF THE ENTROPY 


Ref. in] suggested a thermodynamically motivated 
measure of the gravitational entropy based on the Bel- 
Robinson tensor, 

J . (29) 

A measure of gravitational entropy constructed from this 
tensor was also considered in [531 HD : using an integral 
over conformal time of the super-energy density W de¬ 
fined by 


W = T^,p^n'^n''nPn^ . (30) 

Note that this super-energy density W is observer depen¬ 
dent and non-negative - which is not a problem per se, 
since the entropy is also observer dependent. 

Following Ref. [S] and imposing the following five con¬ 
ditions for the entropy: non-negative, vanishing only if 
Cfivpa = 0, it should measure the local anisotropy of 
the free gravitational field, reproduce the Bekenstein- 
Hawking entropy of a black hole and increase monotoni- 
cally as structure forms, one can define a thermodynami¬ 
cally motivated measure of the gravitational entropy. For 
the case of a perturbed FL spacetime with Euclidean spa¬ 
tial sections, this takes the form [^ 

S'a^D = ^ (“Vt) 


with A a constant, and where we integrate over a comov¬ 
ing volume Vt> ■ Indeed, this differs from the definition 0 
of the relative entropy between two spacetimes. In par¬ 
ticular, it seems to depend on the whole history of the 


spacetime, being present on the left hand side of Eq.(3I) 


a time derivative. But, since the entropy of a spacetime 
with vanishing Weyl tensor is zero, the entropy of a FL 
spacetime vanishes, so that it can also be considered as 
the relative entropy with respect to the background FL 
spacetime. 

We now want to compare this result with our previous 
result for the case of a freely falling observer in a back¬ 
ground spacetime plus first order perturbations. In this 
particular case, one easily concludes that 

W = (32) 


namely the super energy density is equal (up to a con¬ 
stant) to the Weyl scalar, having the part that comes 
from the dual of the Weyl tensor zero contribution. As a 
consequence the gravitational entropy of Ref. [9] reduces 
to 


S'a^o = I 


PI' 


% Jd dr] 




192 


( 33 ) 


We can now compare the two definitions of entropy of 


Eqs. (26) and (33), obtained in Sec. Ill A and Sec. IIIB 


to describe the large-scale structure of the Universe, in a 
ACDM model. 

In the description adopted here, the perturbations are 
stochastic fields, usually with initially Gaussian statis¬ 
tics. It follows that, for example, the spatially averaged 
quantities are also stochastic fields. If W is a function 
of the perturbations and (A) its average on a given do¬ 
main then, from a theoretical point of view, we only have 
access to the distribution of (A), that is to (A), which 
is the ensemble average of (A). Hence, one has to per¬ 
form the ensemble average of the quantities defined in 

Secs. HUD and null 

To compare the definitions of entropy, we have now to 
choose an averaging procedure. In general, one has two 
possibilities. The first is to average over a volume em¬ 
bedded in a spatial hypersurface, for example the one of 
constant redshift z, while the second is to average over a 
two-sphere defined as the intersection of our past or fu¬ 
ture light-cone with this spatial hypersurface. These two 
averaging procedures turns out to be equivalent for terms 
of the kind ((/(‘h))^), with /($) a linear function of the 
gravitational potential. This is due to the fact that in this 
case we average a quantity already of the second order in 
perturbation theory (and the averages can be performed 
only at the background level, if one stops at second or¬ 
der), for a case in which the shape of the domain of in¬ 
tegration is not important (see Eq. (5.3) of Ref. [25] for 
the light-cone averaging case and Ref. [TS] for the spatial 
averaging case). On the other hand, for terms of the kind 
(/($))2 the average depends on the shape of the domain 
of integration and the two procedures give different re¬ 
sults (see Eq. (5.4) of Ref. [l^ for the light-cone averaging 
case and Ref. [53] for the spatial averaging case). There¬ 
fore, it is important to specify which prescription has to 
be used to obtain a physically meaningful result. As a 
guideline we consider the fact that our entropy should 
describe the entropy of the large-scale structure of the 
Universe; that is, it should characterize the Universe as 
a whole and should be averaged over an extended region. 
One can then immediately exclude the light-cone averag¬ 
ing prescription for several reasons. First, the light-cone 
averaging prescription corresponds to an average over a 
two-sphere with dimensions dependent from the value of 
redshift considered, after we fix the observer. For ex¬ 
ample, for a redshift equal to z = 0.01 the region of 
integration (the two-sphere) would be all inside our local 
universe, where non-linearities became extremely impor¬ 
tant. Even worse, the region of integration goes to zero 
at redshift equal to zero. Furthermore, considering dif¬ 
ferent times of observation for a given observer (i.e. dif¬ 
ferent light-cones) the entropy at a given moment would 
be evaluated using different regions of integration. 

As a consequence, it is easy to understand why the 
light-cone averaging prescription cannot be applied in 
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this context. The right choice for the case in consid¬ 
eration is then the spatial averaging prescription. This 
should not surprise us because the entropy is not an ob¬ 
servable in the standard way 

For both definitions of the entropy, obtained in 
Se cs. |III A| and |IIIB[ the key quantity to be evaluated 
is Now, let us consider a peculiar- 

gravitational potential described by a linear power spec¬ 
trum of fluctuations, with the transfer function given in 
Ref. [29] (without the contribution from the baryons) and 
the cosmological parameters from PLANCK [50]. Con¬ 
sidering the Fourier expansion of the first-order gravita¬ 
tional potential we obtain, in Fourier space (and inde¬ 
pendently from the integration domain P; see Ref. [18]) 

_ 1 f fj Jc 

^ y k ), (34) 


with k) = P the power spectrum of the grav¬ 

itational potential. One can then easily see that the in¬ 
tegration in Eq. (34) has an ultraviolet divergence if one 
takes the linear power spectrum defined in Ref. [50]. We 
shall thus consider this as an effective description and 
hereafter we will assume a cut-off fcuv = O.l/iMpc”^ to 
stay within the linear regime and regularize the ultravi¬ 
olet divergence. 

To assume a constant cut-off kuv is enough for our 
purpose. In fact, despite in general the non-linearity 
scale evolve in time, we are mainly interested in a ACDM 
model, starting from when the cosmological constant be¬ 
gins to dominate. For this case the non-linearity scale 
stays nearly constant and our assumption is justified. As 
we will see in the following, this is the case where the 
different entropy proposals here considered behaviour in 
a different way one with respect to the other. 

While the result in Eq. (341 is enough to evaluate the 
gravitational entropy of Eq. (33), for the relative infor¬ 
mation entropy of Eq. (26 1 we have to evaluate also the 
terms present in the second line. In Fourier space the 


terms ((V^4>)^)x) and {didj^d'‘d^^)'D cancel each other, 
while the third term (V^d))!, gives a result dependent 
from the window function used [26] . If we consider a 
top-hat window function of radius R to smooth the field, 
i.e. we have 


W«(|x|) = 0(R-|x|), 

with 0 the Heavyside function, we then have (see, for 
example. Ref. [I 8 ] l 

(V2$)|, = J ^ (35) 


^ See Ref. |27| for the case of cosmological observables, such as the 
second order luminosity distance/redshift relation [28| . where the 
application of the light-cone averaging prescription is necessary. 


with ji{x) = l/a;^(sin(x) — a;cos(a:)) a spherical Bessel 
function. As can be shown, comparing Eq. (351 with re¬ 
spect to Eq. (341, the contribution of (V24>)|, to the total 
value of the relative information entropy is negligible as 
soon as we consider a window function with, for example, 
a radius R at least one order of magnitude larger than 
the cut-off scale kuv- This is a well motived physical 
choice for the case under consideration because we want 
to stay inside the linear regime and the ultraviolet cut¬ 
off chosen is, indeed, the present threshold of the linear 
regime. Therefore, we will neglect this term hereafter in 
the evaluation of the information entropy of Eq. (26). 


A. Relative Information Entropy 


Let us begin with some general considerations. 
Considering Eq. (26), we have that Syu^v/Vv and 
{Cp,u\pC^'^^^)T> are proportional only if the terms in the 
second line of Eq. (26) give negligible contribution. As 


showed in the previous section, this is the case for the 
case of a Universe described by a perturbed FL space- 
time when we average over a large window function. 


As a consequence, using Eq. (26), we can conclude that 


^ c. Mpi “ 
Vp 






(36) 


is a good approximation. Then, using the result of 
Eq. (341, we determine the behaviour of the relative in¬ 


formation entropy per unit comoving volume Vp. It is 
depicted on Fig. for a CDM model (left panel) and a 
ACDM model (right panel). As shown by these two fig¬ 
ures, the expression of Eq. (36) is monotonically increas¬ 


ing with the time only for a CDM model, i.e. as long as 
the cosmological constant is vanishing. The fact that the 
relative information entropy is not a valid definition of 
entropy for a ACDM universe can be understood by the 
fact that in the limit for which the proper time goes to in¬ 
finity (equivalent to z —^ — 1 ) the contribution of the cos¬ 
mological constant is more and more dominant. The Uni¬ 
verse, both at the background level and at the perturbed 
level, is attracted toward a de Sitter spacetime. There¬ 
fore, the relative information entropy between these two 
spacetimes stops growing and then decreases asymptoti¬ 
cally to zero because all scalar perturbations are washed 
out. The turn-over corresponds roughly at the time when 
the cosmological constant starts dominating the cosmic 
expansion. 


To conclude, it is easy to see from Eqs. ( |36[ ) and ( |34[ ), 
that the relative information entropy evolves as ~ 0 ^( 77 ) 
for the CDM case, as already stated in Ref. [13] . 

Let us stress that the shape of the relative information 
entropy in Fig. [^ and of the gravitational entropy in the 
figures in the follow, are independent from the value of 
the ultraviolet cut-off used. This is a consequence of the 
fact that we use a power-spectrum with a linear transfer 
function (as given in Ref. [53]) for which the k and time 
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dependence are factorized. If we had considered a non¬ 
linear transfer function (like the ones in Refs. (STJ |32]) 
then the k and time dependence cannot be factorized 


anymore, and the shape of our figures would be slightly 
dependent from the value of the choice of the ultraviolet 
cut-off. 


COM 


ACDM 



1/(1 + z ) 



FIG. 1: Evolution of the relative information entropy per unit comoving volume, obtained from Eqs. (361, for the CDM model 
(left panel) and a ACDM model (right panel), setting SttG = 1 and assnming an ultraviolet cut-off fcuv = O.l/iMpc”^, as a 
function of 1/(1 -I- 2 ). 


B. Gravitational Entropy 

Let us now perform a similar analysis for the grav¬ 
itational entropy defined in Ref. |9], and described in 
Sec. irTml 

For a perturbed FL universe, at linear order, the ex¬ 


pression (331 does not take into consideration the fun¬ 


damental fact that the perturbations of our Universe are 
stochastic fields. To that purpose, we need to perform 
both a spatial and an ensemble average. Since the Weyl 
scalar is second order in perturbation, we 

need to first average it before we insert it into the en¬ 
tropy expression. This corresponds to modify Eq. (33) 


as 


S\ 


G,V 


a d 




192 


(37) 


where A is a numerical constant that, similarly to what 
done in Ref. [9] to recover the Bekenstein-Hawking en¬ 
tropy for a stationary black hole, we fix equal to one 
hereafter. Using the general solution of Eq. (34), we 


easily obtain the behaviour of the volume entropy for a 
CDM and a ACDM model. They are depicted on Fig. 
for several values of the cosmological constant. In all 
cases the gravitational entropy (37) is monotonically in¬ 


creasing with the time. In particular, in a CDM model, 
the entropy goes to infinity, but with a different time be¬ 
haviour as compared to the relative entropy (361, namely 
it behaves as ~ a^/^(? 7 ). 


In the presence of a cosmological constant the entropy 
tends to a constant in the limit for which the proper 
time goes to infinity. From Fig. one can see how the 
entropy asymptotically reaches its constant value, that 
indeed depends on the value of the cosmological constant. 
The maximum entropy decreases with the cosmological 
constant. The dependence of the maximum asymptotic 
value of the gravitational entropy is plotted as a function 
of the cosmological constant in Fig. [^and can be shown 
to roughly behave as I/Uaq. 


Let us try to explain the physical reasons of such a 
behavior for the entropy in the presence of a cosmolog¬ 
ical constant. The entropy (37) is the entropy associ¬ 
ated to the formation of the large-scale structure, and 
is defined as an integrated effect over the cosmological 
history. The entropy increases when more structures are 
formed. When the Universe approaches the de Sitter 
phase, the growth of structures freezes and so does the 
entropy. Therefore, the gravitational entropy encodes the 
fact that the Universe is asymptotically de Sitter, but it 
only includes the entropy associated with the formation 
of structures. 


To better show this last point, we finally plot the 
evolution of the derivative of Sg,vIVv with respect to 
y = 1/(1 -I- 2 ) in Fig. The curve has a turning point 
when the cosmological constant starts to dominate the 
expansion of the Universe and then goes to zero. As a 
consequence, the integrated effect stops. 





















ACDM 


ACDM 




FIG. 2: Value of the gravitational entropy ( |37| l per unit co¬ 
moving volume for a ACDM model and for different values 
of Hao, setting 8nG — 1 and assuming an ultraviolet cut-off 
fcuv = O.l/iMpc”'^, as a function of 1/(1 -I- z). From top to 
bottom, Hao = 0,0.05,0.35,0.68 (standard ACDM model), 
and 0.95. 


FIG. 4: Value of the gravitational entropy derivative with 
respect to y = 1/(1 -|- z) per unit comoving volume, obtained 
using Eq.(37l, for a ACDM universe, setting 8ttG = 1 and 
assuming an ultraviolet cut-off fcuv = O.l/iMpc”^, in function 
of y. 



FIG. 3: Asymptotic value of the gravitational entropy ( |37[ ) 
per unit comoving volume for a ACDM universe as a function 
of the cosmological constant, setting 8nG = 1 and assuming 
an ultraviolet cut-off fcuv = O.l/iMpc”^. 


V. DISCUSSION 


In this manuscript we have computed the entropy of 
the large-scale structure of the Universe starting from the 
proposed definitions of Refs. [9] and m- Our results 
are valid in the regime for which the Universe can be 
described as a perturbed FL spacetime. This is valid at 
late times on scales ~ 10/i“^Mpc or more, where the 
linear regime holds. 

We have argued that the entropy should be evaluated 
using a spatial averaging procedure. This is because the 
entropy is not an observable in the standard way. It 
describes the large-scale structure of the Universe and 
should arise from an average made over an extended re¬ 
gion. 

Starting from the definition of Ref. [13] we obtain the 
relative information entropy (36). As shown in Fig. 
the expression (36) is monotonically increasing with time 


for a CDM model, while it does not satisfy the Penrose 
conjecture 130 for a ACDM model. Therefore, this def¬ 
inition does not seem to be a valid definition of entropy 
for a ACDM universe. In fact, in the limit for which the 
proper time t —> -|-oo, the importance of the cosmological 
constant increases and both FL and perturbed FL uni¬ 
verses can be well approximated by a de Sitter universe. 
Therefore, the relative information entropy stops grow¬ 
ing at some stage and then asymptotically decreases to 
zero. 


On the other hand, starting from the definition [S] we 
obtain the gravitational entropy (37). In this case, the 
entropy always grows with time, both for a CDM and 
a ACDM model, hence satisfying the Penrose conjec¬ 
ture |3IH]- While it tends to an infinite value for a CDM 
model, it saturates to a constant value in the presence of 
a cosmological constant. The entropy (37) is the entropy 
associated to the formation of the large scale structure 
of the Universe. 


The two proposals of Ref. 0 and Ref. m were al¬ 
ready compared in the literature in the case of a LTB 
dust model in Ref. m- It was shown how in both cases 
the entropy of a LTB dust model grows with time satis¬ 
fying the Penrose conjecture. Such models could be used 
only to describe the local universe and, therefore apply 
in this restrict regime. Our results instead apply on large 
scale, where the linear and mildly non-linear regimes are 
valid and the Universe can be describe by a perturbed 
FL spacetime. The entropy obtained here can be seen as 
the entropy of the large-scale structure of the Universe 
and is complementary to the one obtained in Ref. HD. 

To conclude, in jSSj the proposal of Ref. [S| was stud¬ 
ied in a non-perturbative context, with local voids of 
50 — 100 Mpc described by spherical under-dense regions 
with negative spatial curvature and dynamics determined 
by LTB dust models. The results obtained in [35] for 
the gravitational entropy are in good quantitatively and 
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qualitatively agreement with the result presented here, 
when the LTB evolution is in its linear regime. In partic¬ 
ular, the gravitational entropy has an asymptotic behav¬ 
ior similar to the case of the gravitational entropy of the 
large-scale structure for a ACDM model. Therefore, the 
results of [33j provides an important local physics con¬ 
nection to the large scale linear regime described here. 
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It can first be shown that there is only 6 gauge invari¬ 
ant degrees of freedom usually defined as 

4/ = V^+^AF+^(/3 + F'), (Bl) 

$ ^ a-|(/3 + F0-i(/3 + F7, (B2) 

¥ = +Ib^ , (B3) 

where the prime denotes the derivative with respect to 
conformal time and % = a'/a. 

Considering a matter sector described by a perfect fluid 
with stress-energy tensor 

(p -f P)u^u^ + , (B4) 

where the density and pressure can be split as p{r], x) = 
p(o)(jy) _|_ p(P(^r],x) and P{r],x) = P^^\r]) + P^^'>{ri,x), 
and the velocity of the comoving observers is decomposed 
as + 5u^ with = —1. It follows that 

—a,v'‘) , = a{—l —a,Vi — l/2Bi) (B5) 


Appendix A: Dynamics of the background spacetime and we decompose v, into scalar and a vector component 

according to 


We assume that the background spacetime is well- 
described by a spatially Euclidean Friedmann-Lemaitre 
universe with the late time dynamics dictated only by 
pressureless matter and cosmological constant with den¬ 
sity parameters 


Vi = diV + Vi 


(B6) 


The scalar shear, given in Eq. (10), is by construction 
a second order quantity so that 


_ STrCpmo 

3iJo 


f^AO 


A 


(Al) 


that satisfy flmo + ^ao = 1- The Friedmann equation 
then takes the usual form 


(cr2)(0) ^ (^2^(1) ^ (-gy) 

Thus, at the lowest order and for a general metric, we 
have 


H\z) 

HI 


— f^mo(l + + ^^AO ■ 


(A2) 


and, using the proper time t, its solution is given by 

a{t) oc sinh^/^ (^^V^AoHot^ ■ (A3) 

The normalization to the Hubble constant today, Hq, 
implies that 


sinh Q-\/OAoi?o^o^ = 


n 


1/2 

AO 


_ 3/2 

so that the redshift is given by 

l + z= ^0 

sinh^/^ {^^/n AO Hot) 

Appendix B: Linear perturbation theory 


(A4) 


(A5) 


The evolution of the degrees of freedom of the met¬ 
ric (13) can be found in many textbooks, e.g. Ref. [T] . 


(a2)(2) 


1 

2a^SB 



6S^,j6SH - i(v25S') 


2 


1 

2a‘^S' 






1 




_ B h 

2o? ^'^ 


+ (B8) 

where h,ij = \DijE + d(iXj) + \hij, we use the notation 
Xi = diX for any field X, and 6S is the first order 
perturbation of the scalar <5'(x, t) defining the space-time 
foliation. 


It is clear that first order perturbation theory is suffi¬ 
cient to obtain the general expression for the shear up to 
second order, since second order perturbations will con¬ 
tribute only to third or fourth order to tr^. 
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